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Abstract: Plants are an important component in many natural ecosystems. They influence soil properties, 
especially in arid ecosystems. The selection of plant species based on their adaptations to site conditions is 
essential for rehabilitation of degraded sites and other construction sites such as check-dams. Other 
factors to be considered in species selection include their effects on soil properties and their abilities to 
meet other management objectives. The purpose of this study was to assess the effects of native (Populus 
euphratica Oliv. and Tamarix ramosissima Ledeb.) and introduced (Eucalyptus camaldulensis Dehnh. and Prosopis 
juliflora (Swartz) DC.) woody species on soil properties and carbon sequestration (CS) in an arid region of 
Iran. Soil sampling was collected at three soil depths (0—10, 10—20 and 20—30 cm) at the sites located 
under each woody species canopy and in an open area in 2017. Soil physical-chemical property was 
analyzed in the laboratory. The presence of a woody species changed soil characteristics and soil CS, 
compared with the open area. For example, the presence of a woody species caused a decrease in soil bulk 
density, of which the lowest value was observed under E. camaldulensis (1.38 g/cm?) compared with the 
open atea (1.59 g/cm?). Also, all woody species significantly increased the contents of soil organic matter 
and total nitrogen, and introduced species had more significant effect than native species. The results 
showed that CS significantly increased under the canopy of all woody species in a decreasing order of P. 
euphratica (9.08 t/hm?)>E. camaldulensis (8.37 t/hm?)>P. juliflora (5.20 t/hm?) T. ramosissima (2.93 
t/hm?)>open area (1.33 t/hm?), thus demonstrating the positive effect of a woody species on CS. 
Although the plantation of non-native species had some positive effects on soil properties, we 
recommend increasing species diversity in plantations of native and introduced woody species to provide 
more diversity for the increased ecosystem services, resilience, health and long-term productivity. 
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1 Introduction 


Arid regions occur where evaporation is higher than the rainfall rate (Gaur et al., 2018). More 
than a third of the Earth's population (2.5x10? people) live in arid regions, which cover about 
5.36x10° km? of the Earth (41% of the surface area) (Mortimor et al., 2009; UNEMG, 2011). In 
recent decades, these arid regions have experienced accelerated development, along with changes 
in land use and climate (Song et al., 2018). As a consequence, arid regions are increasingly being 
degraded, and the annual index of desertification is 5.8x10° hm?/a (Lal, 2001). In desert 
ecosystems, soil carbon content decreases and the emission of carbon as CO» into the atmosphere 
accelerates (Lal, 2004; Romm, 2011), thus contributing to the global warming. In addition, 
desertification often leads to soil erosion and salinization, losses of vegetation cover and biomass, 
and decreases in soil productivity and quality (Zhao et al., 2006; Ma et al., 2017). 

The selection of appropriate woody species in reclamation and rehabilitation projects in arid 
regions is important because woody species have a great influence on the degree of success in the 
recovery of ecosystems, and they increase productivity in goods and ecological services for the 
people. Some researchers have been done in revegetation of arid regions such as Kergoat et al. 
(2018) and Verón et al. (2018) who assessed the quantity and type of vegetation cover for 
preventing desertification, soil erosion and ecosystem deterioration in arid areas. Nyssen et al. 
(2009) showed that establishing eucalypt woodlands had a positive influence on the environment, 
especially on ameliorating soil conditions. Tesfay et al. (2015) compared the ability of different 
tree species to improve soil conditions and produce wood fuels. And Reubens et al. (2011) used 
Decision Support Systems to identify appropriate tree species for rehabilitation of arid regions 
and prioritize woody plants for revegetation and reclamation based on their adaptations to the 
environment. However, selection of tree species is complicated as it is influenced by reclamation 
goals, existing vegetation, degree of site degradation and ecosystem processes (Reubens et al., 
2011). To date, few studies have been done to evaluate appropriate plant species for reclamation 
projects in certain arid regions, including western Iran. 

Iran is centrally located in the arid and semi-arid regions of the Earth, with more than 6096 of 
its land area (about 105x109 hm?) being classed as arid and semi-arid (Modarres and da Silva, 
2007). In Iran, woody species plantations are commonly used to amend the soil, combat 
desertification, revegetate large construction sites and mitigate dust storms. Choosing suitable 
woody species that are adapted to and can tolerate hot and dry conditions are needed. For example, 
eucalyptus (Eucalyptus camaldulensis Dehnh.), Euphrates poplar (Populus euphratica Oliv.), 
saltcedar (Tamarix ramosissima Ledeb.) and mesquite (Prosopis juliflora (Swartz) DC.) are 
among the most commonly used woody species in these arid and semi-arid regions in Iran, in 
which P. euphratica and T. ramosissima are native species, and E. camaldulensis and P. juliflora 
are introduced species. Although the field performance of these woody species, i.e., survival and 
growth is good, the adapted strategy about how different origins of these tree species modify the 
environment and the soil is lack. Therefore, the objective of this study was to assess the effects of 
introduced and native woody species on soil properties and carbon sequestration (CS) in arid 
regions. Our results may be useful to managers and decision makers who need to select woody 
species and optimize their efforts in rehabilitating degraded arid ecosystems. 


2 Materials and methods 


24 Study area 


The study area in the Reza Abad Park (46?13'-46?14'N, 33?09'-33?12'E) is 700 hm? and is 
located in the Mehran County of Ilam Province in the west of Iran (Fig. 1). The mean annual 
precipitation for the area is 209 mm, which mostly occurs in the autumn and winter seasons. The 
annual mean temperature is 24.5?C. According to the climatic curve (precipitation and 
temperature), ten months of the year are dry (Fig. 2) with a warm and dry climate. The elevation 
of the study area ranges from 203 to 253 m above sea level, and the land is relatively flat with less 
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than one percent slope. The soil is calcareous with more than 3096 gravel on the surface, low in 
organic matter and shallow with a deep ground water table (Zamani et al., 2018). The dominant 
soil texture is sandy loamy. The geology of the study area is of the Gachsaran formation, which is 
characterized by alternating layers of anhydrite gypsum, marl and claystone. 


46?12'15"E 46?14'15"E 


33?09'40"N — 33?11'10"N 


46?12'I5"E . 46?14'15"E 


Fig. 1 Location of the study area in Ilam Province, Iran 
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Fig.2 Mean monthly precipitation and monthly mean temperature in the study area 
2.2 Woody species 


Study sites included a natural forest stand of native T. ramosissima and P. euphratica species 
along the Konjancham River, and 26-year-old plantations of introduced E. camaldulensis and P. 
juliflora species in the Reza Abad Park (Table 1). The common understory vegetation was 
comprised of Stipa capensis, Plantago ovata and Anchusa strigose. 


Table1 Stand characteristics of woody species in the study area 


Forest type Woody species Age (a) Height (m) Canopy cover (m°) DBH (cm) 
Introduced P. juliflora 25-27 4.0+0.4 12.641.3 15.4+0.8 
species E. camaldulensis 25-27 9.7+0.7 5,141.4 13.5+0.5 
T. ramosissima - 2.7+0.2 3.50.2 3.20.3 


Native species 
P. euphratica - 3.2+0.4 7.8t1.7 10.5+0.2 


Note: - means no data. DBH, diameter at breast height. 


2.3 Soil sampling design 


At each study site, we randomly established eight 20 mx20 m plots. In September 2017, soil 
samples from three different depths (0—10, 10-20 and 20-30 cm) were collected from each plot 
under the woody canopy and in the open area (control). Each sample was a composite of four 
sub-samples taken around one species at the middle of the plot in the forest and plantations. In 
addition, one intact soil core (5-8 cm height and 8 cm diameter) was taken at each sampling plot 
for bulk density (BD) determination. 


2.4 Soil physical and chemical analyses 


Soil samples were treated with a 2-mm sieve to removed roots and coarse fragments and stored in 
plastic bags at room temperature before analyses. Soil BD was measured by common core method 
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(Black and Hartge, 1986). Soil porosity (n) was calculated from BD and particle density (Weil 
and Brady, 2015). Saturation point (SP) was determined gravimetrically. Field capacity (FC) at 33 
kPa and permanent wilting point (PWP) at 1500 kPa were determined using a pressure plate 
apparatus (Klute, 1986). Soil pH was measured with a glass electrode in a soil suspension (1:1 
w:v). Soil electrical conductivity (EC) was measured with an EC meter. The Walky and Black wet 
oxidation method was used to determine soil organic carbon (SOC) content (Nelson and Sommers, 
1986), which was also used to calculate organic matter (OM) by multiplying 1.724 times the 
organic carbon value (Weil and Brady, 2015). The Kjeldahl procedure was used to determine total 
nitrogen (TN) content (Bremner and Mulvaney, 1982). Cation exchange capacity (CEC) was 
measured through sodium acetate replacement (pH-8.5) (Summer and Miller, 1996). Soil 
available phosphorus (AP) content was determined by NaHCO; extraction (Olsen and Sommers, 
1982). Exchangeable cations of Ca and Mg from three soil depths were extracted using a 0.1-M 
BaCl2 (Hendershot and Duquette, 1986). Equations 1 and 2 were used for estimating total 
dissolved salts (TDS) (Chang et al., 1983) and CS (Qin et al., 2016), respectively: 


TDS=765.1 EC!097, (1) 
CS-SOCxBDxsoil depth. (2) 
2.5 Statistical analyses 


Prior to statistical analyses, we analyzed the normality of means and homogeneity of variances 
using Kolmogorov-Smirnov and Levene tests, respectively. For the statistical analyses, we used a 
generalized linear mixed effect model (GLMM) with a Poisson error distribution and a log-link 
function where, species type, depth and their interaction were fixed factors and individual species 
was a random factor in the model. GLMM is an extension of the class of generalized linear model 
that adds a random effect to the linear predictor (McCulloch and Neuhaus, 2005). 

GLMM analysis was separately done for all soil physical, chemical and carbon sequestration 
response variables. This analysis was performed using the 'Ime4' package (Bates et al., 2015). 
Stepwise regression was used to determine the main factors affecting CS. Before regression 
analysis was done, the collinearity among independent variables (soil physical and chemical 
properties) was evaluated using Pearson's correlation coefficients. 

In addition, we removed BD and SOC from the final model because they were used to calculate 
CS. The means comparison models (i.e., ANOVA (analysis of variance), GLMM or mixed models) 
only evaluate a single soil property, and lack the ability to provide for the general compression of 
multivariate overall soil structure among treatments. 

Therefore, the nonmetric multidimensional scaling (NMDS) test using the Bray-Curtis 
similarity matrices was applied to the overall structural changes in soil physical and chemical 
properties. After conducting NMDS, we used ANOVA to determine the significant differences 
among groups of the NMDS axis (Yang et al., 2019). The NMDS analysis was performed using 
the 'vegan' package (Oksanen et al., 2018). All analyses were performed with the last version of R 
statistical software (R Core Team, 2018). 


3 Results 
3.1 Soil physical properties 


Soil physical properties were significantly affected by woody species, and soil depth had a 
significant effect on some properties such as BD, n and FC (P«0.05; Table 2). The interaction 
between woody species and sampling depth did not have a significant effect on any soil physical 
properties. In other words, variations in soil physical properties under different woody species 
were independent from sampling depth. 

Woody species caused a decrease in BD based on means comparison tests on soil properties. 
BD was significantly lower (P=0.002) under E. camaldulensis (1.38 g/cm?) than in the open area 
(1.59 g/cm5; Fig. 3). BD significantly increased with increasing soil depth from the first soil depth 
(0—10 cm) to the third soil depth (20-30 cm), while in the open area there was no significant 
difference in BD among the three soil depths (Fig. 3). 
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Table 2 Results from generalized linear mixed-effects (Chi-square values and P-values) model with a Poisson 
family distribution for the effects of species, depth and their interaction on soil properties 


Species Depth SpeciesxDepth 
Soil property Variable 

x P-value x P-value x P-value 
BD 16.70 0.0020 54.470 «0.0001 6.12 0.1910 
n 16.42 0.0020 54.530 «0.0001 6.07 0.1930 
Physical SP 294.42 «0.0001 3.380 0.0657 1.41 0.8420 
FC 126.08 «0.0001 6.240 0.0120 4.33 0.3620 
PWP 74.15 «0.0001 2.770 0.0950 1.96 0.7460 
pH 515.88 «0.0001 6.190 0.0120 20.50 «0.0001 
EC 200.20 «0.0001 17.020 «0.0001 12.65 0.0130 
OM 510.44 «0.0001 14.560 «0.0001 4.33 0.3620 
TN 217.21 «0.0001 127.570 «0.0001 42.44 «0.0001 
i AP 6.48 0.1660 1.290 0.2560 6.90 0.1410 
Ca 1365.01 <0.0001 1069.000 0.0010 9.80 0.0430 
Mg 1232.08 <0.0001 6.270 0.0120 23.79 <0.0001 
CEC 519.64 <0.0001 0.003 0.9520 8.50 0.0740 
TDS 361.12 <0.0001 30.070 <0.0001 11.69 0.0190 
CS 612.13 <0.0001 0.770 0.3800 2.64 0.6190 


Note: BD, bulk density; n, porosity; SP, saturation point; FC, field capacity; PWP, permanent wilting point; EC, electrical conductivity; 
OM, organic matter; TN, total nitrogen; AP, available phosphorus; Ca, exchangeable calcium; Mg, exchangeable magnesium; CEC, 
cation exchangeable capacity; TDS, total dissolved salts; CS, carbon sequestration. 


The value of n increased under all woody species, but it was significantly higher under E. 
camaldulensis (20%) than in the open area. The highest value of n was observed at the 0-10 cm 
depth (5396), whereas the lowest value of n occurred at the 20-30 cm depth (38%). There was no 
significant variation in the value of n with soil depths in the open area. SP was significantly 
higher under P. euphratica and T. ramosissima species (3396 and 1696, respectively) than in the 
open area. In contrast, it was not significantly different between the open area and under E. 
camaldulensis or P. juliflora species. SP did not vary significantly with soil depths (Fig. 3). FC 
was significantly higher under P. euphratica and T. ramosissima species (3896 and 2996, 
respectively), but not under E. camaldulensis or P. juliflora species than in the open area. 
Sampling depth had a significant effect on FC only under P. juliflora species, while under the other 
woody species, there was no significant difference among different soil depths. Finally, the 
canopy of all woody species caused a significant increase in PWP, and the highest values of PWP 
were observed in the soils under the canopies of P. juliflora (7.4396) and E. camaldulensis 
(7.0896). 


3.2 Soil chemical properties 


The results of GLMM analysis showed that all the soil chemical properties except AP (P=0.166) 
were significantly affected by woody species. Also, soil chemical parameters were significantly 
affected by soil depths except for CEC (P«0.001) and CS (Table 2). The interaction between plant 
species and sampling depth was also significant for pH, EC, TN, Ca, Mg and TDS parameters 
(P«0.05), indicating the dependence of the changes of these factors on the soil sampling depth. 
Soil pH significantly increased under all woody species except for P. euphratica, compared with 
the open area (Fig. 4). Soil depth also had a significant but variable effect on pH depending on 
woody species. For example, soil pH significantly increased with increasing soil depths under P. 
juliflora species (0—10 cm, pH=7.49; 10-20 cm, pH=7.50; and 20-30 cm, pH=7.57), but under T. 
ramosissima species, it decreased with increasing soil depths (0-10 cm, pH=7.34; 10-20 cm, 
pH=7.35; and 20-30 cm, pH-7.27). 

The relationship between EC and plant species varied such that EC values of E. camaldulensis 
and T. ramosissima increased up to 1796 and 2796, respectively, but it decreased significantly up to 
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5896 for P. juliflora, while there was no significant difference in EC between P. euphratica and the 
open area. EC increased with increasing soil depth under E. camaldulensis and P. euphratica 
species (Fig. 4). 

OM content significantly increased under P. euphratica followed by E. camaldulensis, P. 
juliflora and T. ramosissima (8396, 8296, 7096 and 5496, respectively) compared with the open 
area. OM content was higher at the soil surface layer than at the other soil depths. CEC content 
significantly differed (P«0.001) among woody species, and it was significantly higher under any 
woody species compared with the open area (6.71 cmol/kg) with the exception of T. ramosissima 
(7.46 cmol/kg). CEC did not vary significantly at different soil depths. 
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Fig. 3 Soil physical properties (a—e) under native (Tamarix ramosissima and Populus euphratica) and introduced 
(Eucalyptus camaldulensis and Prosopis juliflora) woody species and in the open area (control). Different 


lowercase letters indicate significant difference among different soil depths at P<0.05 level, and different uppercase 
letters indicate significant difference between woody species and the open area at P<0.05 level. 


TN was significantly higher under all woody species than in the open area (0.001%), except for 
T. ramosissima (0.007%). TN was also significantly higher in the upper soil depth (0-10 cm) 
under all woody species compared with the other soil depths. Both T. ramosissima and P. 
euphratica significantly increased Mg (69% and 8%, respectively) and Ca contents (55% and 54%, 
respectively) of the soil compared with that of in the open area. Soil CS under all woody species 
significantly decreased in the following order P. euphratica (9.08 t/hm?)* E. camaldulensis (8.37 
t/hm?)>P. juliflora (5.2 t/hm?)>T. ramosissima (2.93 t/hm?)»open area (1.33 t/hm?), indicating the 
positive effect of a woody species on CS. There was no significant difference in CS among 
different soil depths (Fig. 4). 
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Fig. 4 Soil chemical properties (a-h) and carbon sequestration (CS, i) under native (Tamarix ramosissima and 
Populus euphratica) and introduced (Eucalyptus camaldulensis and Prosopis juliflora) woody species and in the 
open area (control). Different lowercase letters indicate significant difference among different soil depths at 
P<0.05 level, and different capital letters indicate significant differences between woody species and the open 
area at P«0.05 level. EC, electrical conductivity; OM, organic matter; TN, total nitrogen; Ca, exchangeable 
calcium; Mg, exchangeable magnesium; CEC, cation exchangeable capacity; TDS, total dissolved salts. 
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3.3 Factors effecting soil CS 


Individual model predicting CS content of treated woody species showed that soil CS was 
affected by FC under E. camaldulensis species, by Mg and n under P. juliflora species, and by FC 
and SP under P. euphratica species. In contrast, soil CS was not affected by any of the physical 
and chemical variables of the soil under T. ramosissima and in the open area. Finally, regression 
model showed that soil CS was affected by TN, Mg, TDS, CEC and PWP (Table 3). 


Table 3 Results of multiple regressions between soil carbon sequestration and different physical and chemical 
soil properties in each treatment separately and for the combined data 


Treatment Independent variable Final model r R? P-value 
Control" - - - - - 
E. camaldulensis FC —0.419(FC)+13.111 0.604 0.365 0.037 
P. juliflora Mg, n —0.092(n)+9.4369 0.733 0.537 0.007 


T. ramosissima” - - - - - 


P. euphratica FC, SP 0.736(FC)-1.826 0.731 0.535 0.007 
Overall data TN, Mg, TDS, CEC, PWP.  69.149(TN)*2.282(TDS)-0.075(Mg)*0.506 0.918 — 0.842 0.000 


Note: * means that no predictor variables were interred to the model. FC, filed capacity; Mg, exchangeable magnesium; n, porosity; SP, 
saturation point; TN, total nitrogen; TDS, total dissolved salts; CEC, cation exchangeable capacity; PWP, permanent wilting point; -, no data. 


3.4 NMDS test 


The results showed that soil properties can be used to differentiate among the treatments, both in 
the open area and under woody species sites. One way ANOVA result showed that there was a 
significant difference between the treatments in the first axis (F-value=103; P-value=0.000) and 
those in the second axis (F-value=250.6; P-value=0.000). The most important factors 
characterizing the habitat of T. ramosissima were EC, SP, Ca, BD, Mg and FC. In contrast, the 
factors that best defined the habitat of P. juliflora were acidity of the exchange capacity of 
phosphorus and n. The environmental requirements of both P. euphratica and E. camaldulensis 
are similar in terms of OM, TN and soil CS levels and hence, they had the least distance between 
each other. High values of TN, OM and soil CS (mostly in P. euphratica) and PWP (mostly in E. 
camaldulensis) characterized these stands. The control area is located far from the study forest 
and plantations. According to the orientation of the axes, we concluded that the minimum 
amounts of OM, CS, and TN found in the open area supported the importance of the plant species 
in improving soil conditions (Fig. 5). 


4 Discussion 


4.4 Comparing soil physical properties in native and introduced woody species 


Our results demonstrated that woody species on the native or introduced plantation areas generally 
affected the suite of soil physical properties including BD, n, SP, FC and PWP. These plant species 
effects on soil physical properties were most noticeable when compared with those in the open 
area. Our results are consistent with other studies (Frouz et al., 2013; Kooch et al., 2016), 
especially in arid and semi-arid regions (Kalinda et al., 2015; Chen et al., 2016). The lowest BD 
was observed under E. camaldulensis species, which also had a stronger influence on n compared 
with that of in the open area. In general, plant litterfall increases OM input to the soil (Carnol and 
Bazgir, 2013), similar to the contributions from tree roots. Tree roots also create numerous pores 
in the soil increasing n and decreasing BD (Binkley and Fisher, 2012). In addition, a tree canopy 
has a positive effect on the activity of mesofauna and macrofauna of the soil by improving 
moisture in the soil environment and providing underground energy sources (Heydari et al., 2017; 
Zagatto et al., 2019). The activity of soil fauna is an important factor that helps to reduce BD by 
creating pores and displacing OM. 

The effect of E. camaldulensis on decreasing of BD is stronger than those of the other species 
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Fig. 5 First two dimensions nonmetric multidimensional scaling (NMDS) ordinations based on the Bray-Curtis 
dissimilarity matrices of soil physical and chemical properties. BD, bulk density; n, porosity; SP, saturation 
percentage; FC, filed capacity; EC, electrical conductivity; PWP, permanent wilting point; OM, organic matter; 
TN, total nitrogen; AP, available phosphorus; Ca, exchangeable calcium; Mg, exchangeable magnesium; CEC, 
cation exchangeable capacity; TDS, total dissolved salts; CS, carbon sequestration; Ta, T. ramosissima; Po, P. 
euphratica; Eu, E. camaldulensis; Pr, P. juliflora. 


may be due to a more extensive root system and higher levels of OM production (leaves and bark of 
the trunk and branches). In contrast, Montero and Delitti (2017) reported that BD increased under E. 
camaldulensis due to soil compaction. In this study, BD increased and n decreased with soil depths 
under woody species, while there was no significant change with soil depths in the open area. 
Basically, BD increased with soil depths because of the weight of the upper layers and lower OM 
content resulting from more compaction and lower soil porosity (Weil and Brady, 2015). 

The SP and FC under P. euphratica species was the highest among the woody species and the 
open area. This may be due to an increase in soil OM from the understory trees (Hofstede et al., 
2002; Farley et al., 2004). The stability of SP and FC with soil depths may be due to the uniformity 
of soil texture with depths in the study area (Zhao et al., 2011; Weil and Brady, 2015). The results 
showed that AP was not influenced by the introduced or native woody species. Similar results 
were reported by Isichei and Muoghalu (1992) and Son et al. (1992) who found that AP was not 
affected by plantations. Phosphorus has a low mobility in the soil (Havlin et al., 2016), and other 
factors such as climate have more effect on changes in soil P than vegetation. Heydaei et al. (2019) 
reported that soil phosphorus was significantly influenced by the climate, but not by management. 
The lack of differences in phosphorus between plantations and natural forest may be because it is 
not limiting in the study area. 


4.2 Comparing soil chemical properties in native and introduced woody species 


All soil chemical properties including pH, EC, SOM, TN, exchangeable Mg and Ca, CEC and 
TDS were affected by woody species, demonstrating the important influence of plant species on 
soil chemistry. Many studies have reported the effect of the tree canopy and shrub species on soil 
chemical properties (Yang et al., 2011; Waring et al., 2015; Habashi et al., 2019). Soil depth also 
has a significant effect on most soil chemical properties (Weil and Brady, 2015). Bardgett (2005) 
and Binkley and Fisher (2012) reported that root secretion and soil microbial biomass had more 
effect on chemical properties than physical properties (Bardgett, 2005). Movement of nutrient 
elements and water-soluble materials from the topsoil into the subsoil by leaching leads to the 
differences between soil horizons (Lee and Jose, 2005; Buol et al., 2011). 

Among the studied woody species, E. camaldulensis and P. euphratica had more effect on soil 
chemical properties such as OM, TN, CEC and TDS. Soil has a higher OM under woody species 
than in the open area because of inputs from litterfall and roots (Prescott, 2002). Both E. 
camaldulensis and P. euphratica are fast growing trees with large canopies that produce more 
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litterfall than the other species (Singh et al., 1989; Williams and Wardle, 2007). In addition, E. 
camaldulensis has high contents of tannins and aromatics in the leaves, which increases resistance 
to microbial decomposition, thus resulting in OM accumulation in the forest floor (Coleman et al., 
2004; Brennan et al., 2009). Furthermore, OM from annual bark shedding under E. camaldulensis 
results in the greater input that that of under P. euphratica (Cornelissen et al., 2017). Moreover, 
there is a positive correlation between OM and CEC (Helling et al., 1964; Obalum et al., 2017). 
Thus, these woody species might increase soil fertility and productivity through improving the 
nutrient and water holding capacity and availability in arid and semi-arid regions. The increase in 
soil nutrient accumulation of forest plantations in arid areas might create the effect of fertile 
island. And the degree of development of fertile islands under individual plants depends on the 
canopy size, developmental stage and duration of litter fall accumulation (Li et al., 2008). 

Soil pH was the highest and EC was the lowest under P. juliflora species. Consistent with our 
results, Bruckner (2012) found that there is a negative relationship between soil pH and EC 
because soil pH decreases but the hydrogen concentration increases in the soil solution, which 
leads to an increase in EC. In contrast, EC increased under woody species of E. camaldulensis 
and T. ramosissima. Some species such as T. ramosissima can uptake salt from groundwater, 
transport salt to the leaves, and then return salt to the soil through the litterfall, which increases 
soil salinity (Arndt et al., 2004; Stromberg et al., 2009). Exchangeable Ca and Mg under native 
tree species in natural forests, i.e., P. euphratica and T. ramosissima in this study, were higher than 
those of under the introduced woody species grown in plantations, or in the open area. Leaf 
content of cations such as Ca and Mg varies among woody species, and consequently, the amount 
of nutrients returned to the soil differs by species. For example, Meiresonne et al. (2006) reported 
that a higher exchangeable Ca found under poplar species (P. trichocarpa and P. deltoides) than 
under the other species is a result of the accumulation of this element in the foliar litter of the 
Populus species. In this study, the native trees of the natural forests are older than the introduced 
species grown in plantations; therefore, the higher amounts of aboveground and belowground 
biomasses in the natural forest can result in the higher exchangeable Ca and Mg contents in the 
soils (Binkley and Fisher, 2012). 

We observed that woody species increased soil CS, which was higher under E. camaldulensis 
and P. euphratica species than under the other species. Trees sequester atmospheric CO» via 
photosynthesis and store it in shoot and root carbon-based biomasses. Eventually, the carbon 
enters the soil as OM when the tree sheds or loses leaves, roots, twigs, bark and wood. Therefore, 
CS is higher in natural forests and plantations than in the open areas. In this regard, Grünzweig et 
al. (2003) reported the positive effect of a tree canopy on the increased CS in arid areas. The 
higher amount of soil CS under E. camaldulensis and P. euphratica species may be a result of 
high levels of biomass from high rates of canopy and root growth. Quideau et al. (2001) and 
Pérez-Bejarano et al. (2010) reported that the amount of soil CS is closely related to the plant 
species. In addition, other studies have shown that litterfall inputs not only increase soil carbon 
content, but also promote microbial activities, and hence, nutrient cycling (Xuluc-Tolosa et al., 
2003; HagenThorn et al., 2004). 

FC was the best predictor variable of soil CS under E. camaldulensis and P. euphratica species, 
while n had the greatest impact on the amount of soil CS under P. juliflora. Generally, arid 
ecosystems are characterized by a soil moisture restriction (Gaur et al., 2018). It has been shown 
that water availability and water use efficiency play a critical role in soil CS (de Deyn et al., 2008), 
especially in arid ecosystems. In addition, the traits of complementarity and facilitation in terms 
of water storage and use, as well as protection from drought and solar radiation, are most likely to 
enhance soil CS (Schenk and Jackson, 2002). In other words, plant traits that enable opportunistic 
use of precipitation events play a vital role in soil CS in arid ecosystems (de Deyn et al., 2008). 
Therefore, E. camaldulensis and P. euphratica may more effectively use the available water 
(precipitation and underground water) due to these traits such as a higher leaf area and a 
deep-rooting system. Indeed, the dense canopy coverages of E. camaldulensis and P. euphratica 
likely reduce water and carbon losses (Arriaga and Maya, 2007), which results in optimum water 
use and higher soil SC content. This conclusion is supported by the observed soil carbon content, 
which was the highest under P. euphratica and E. camaldulensis species. 
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Finally, our results demonstrated that native and introduced woody species caused different 
changes in soil properties. Similar results have been reported by other researchers (Pei et al., 2016; 
Wartenberg et al., 2017). Therefore, it is important to select suitable woody species that are 
adapted to improve soil properties and meet other management goals in arid and semi-arid 
regions. 


5 Conclusions 


Selecting the suitable woody species is one way to increase carbon accumulation and soil nutrients. 
Fast-growing species can sequester more carbon in their tissues at a higher rate in a short-term 
than slow-growing species. However, slow-growing species are often longer-lived and they can 
store carbon for longer periods. In this study, fast-growing species such as E. camaldulensis, are 
preferred species to increase soil nutrients and to improve soil properties in arid and semi-arid 
regions because they are capable of producing higher amounts of litter fall and OM to the soils. 
Accordingly, to increase soil CS and soil fertility in arid ecosystems, we suggest establishing 
plantations with species such as eucalyptus and restoration of native forests with species such as P. 
euphratica for management actions. In addition, selection of salt tolerant plants (T. ramosissima) 
to reclaim saline soils could be used as an alternative management approach. It seems desirable to 
increase plant species cover and diversity and to consider the more positive effects of non-native 
species on soil ecology, when designing rehabilitation plans in arid and semi-arid regions. A 
mixed plantation of native and non-native, fast- and slow-growing, short- and long-lived species 
may increase biodiversity, optimize CS and soil improvement, increase resistance and resilience 
to disturbances and environmental stresses, and provide a more diverse mix of ecosystem goods 
and services. 
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